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ABSTRACT 

We present Keck/MOSFIRE /Gbaud spectroscopy of the first mass-selected sample of galaxies at z ~ 2.3. 
Targets are selected from the 3D -Hubble Space Telescope Treasury survey. The six detected galaxies have a mean 
[N ii ]/.6584/Htf ratio of 0.27 ± 0.01, with a small standard deviation of 0.05. This mean value is similar to that of 
UV-selected galaxies of the same mass. The mean gas-phase oxygen abundance inferred from the [N n]/Ha ratios 
depends on the calibration method, and ranges from 12+log(0/H) gas = 8.57 for the Pettini & Pagel calibration 
to 12+log(0/H) gas = 8.87 for the Maiolino et al. calibration. Measurements of the stellar oxygen abundance in 
nearby quiescent galaxies with the same number density indicate 12+log(0/H) stars = 8.95, similar to the gas-phase 
abundances of the z ~ 2.3 galaxies if the Maiolino et al. calibration is used. This suggests that these high-redshift 
star forming galaxies may be progenitors of today’s massive early-type galaxies. The main uncertainties are the 
absolute calibration of the gas-phase oxygen abundance and the incompleteness of the z ~ 2.3 sample: the galaxies 
with detected Ha tend to be larger and have higher star formation rates than the galaxies without detected Ha, and 
we may still be missing the most dust-obscured progenitors. 
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1. INTRODUCTION 

Observations of the elemental abundances of galaxies provide 

information on the build-up of metals in the Universe and on the 
importance of winds and feedback (Dave et al. 2012). Most 

studies find that the mass-metallicity relation evolves with 
redshift, such that at fixed stellar mass, galaxies have lower 

metallicity at earlier times (Erb et al. 2006b; Maiolino et al. 

2008; Zahid et al. 2013; though see also Stott et al. 2013). This 
is consistent with expectations from simple models in which gas 
is gradually enriched by (post-)asymptotic giant branch stars and 

supernovae. 

In addition to measuring their gas-phase metallicities, it is 

also possible to measure the stellar metallicities of galaxies 
(Gallazzi et al. 2005; Panter et al. 2008; Conroy et al. 2013, 
and references therein). As the stellar metallicities reflects the 
gas-phase metallicities at the time of star formation, the com- 
bined measurements of stellar and gas-phase metallicities over 
cosmic time puts powerful constraints on galaxy formation 
models (Bresolin et al. 2009; Sommariva et al. 2012). 

In this Letter we take a step in this direction by comparing 
the stellar oxygen abundances of massive galaxies in the local 
Universe to the gas-phase oxygen abundances of their putative 
progenitors at early times. This comparison should be relatively 

straightforward for massive galaxies, as they formed most of 


their stars at redshifts z > 2 (Thomas et al. 2005; Conroy 
et al. 2013). Therefore, there should be a direct correspondence 
between the gas-phase metallicities of massive galaxies at z > 2 
and the stellar metallicities of their descendants at z = 0. 

This project has recently become possible due to the CAN- 
DELS and 3D -Hubble Space Telescope ( HST) datasets, which 
provide mass-limited samples with accurate redshifts, and to the 
advent of the MOSFIRE spectrograph on the Keck Telescope 
(McLean et al. 2012). Furthermore, accurate stellar abundances 
of individual elements have recently been derived from aver- 
aged spectra of Sloan Digital Sky Survey galaxies of different 
masses (Conroy et al. 2013). 

In order to link progenitor galaxies with their descendants, 
we require that they have the same cumulative number density 
(van Dokkum et al. 2010; Patel et al. 2013; Leja et al. 2013). 
This comparison at a constant number density is preferable to 
comparison at constant stellar mass, as it explicitly takes the 
mass evolution of galaxies into account. 

Throughout the Letter, we assume a Chabrier initial mass 
function (IMF; Chabrier 2003) and a ACDM cosmology with 
Ho = 70 km s _1 Mpc” 1 , Q m — 0.3, and = 0.7. 

2. SAMPLE SELECTION 

We select spectroscopic targets in the UKIDSS-UDS field 
from version 2.1 of the 3D-HST survey catalogs (Brammer 
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Figure 1. Stellar mass evolution is inferred at a constant cumulative number density of 2 x 1 0 4 Mpc~ 3 (dashed horizontal line) from the mass functions of Marchesini 
et al. (2009) (colored lines). The right panel shows the mass evolution at this number density (solid circles) and a quadratic parameterization (solid line, Equation (1)). 
The selection area (±0.2 dex) is shaded in light blue, and target galaxies are marked with crosses. 

(A color version of this figure is available in the online journal.) 


Table 1 

Spectroscopic Targets 


ID 

a 

(J2000) 

<5 

(J2000) 

F140W 

72 a 

log(M*) 

(M@) 

SFR 

(M q yr _1 ) 

A b 

(kpc) 

Z3D-HST 

Zspec 

[Nil] /Ha 

2522 

34.42735 

-5.26471 

22.69 

24.34 

10.61 

242 

5.26 

2.30 

2.315 

0 29 +0 03 
u - zy -0.02 

5698 

34.44652 

-5.24892 

23.55 

26.24 

10.71 

141 

2.82 

2.11 

2.127 

0 3 1 +0 03 
u - :)1 -0.03 

8461 

34.42749 

-5.23648 

23.66 

26.24 

10.78 

142 

3.31 

2.19 

2.299 

n 9Q+0.03 
0.02 

10746 

34.39999 

-5.22679 

22.39 

24.09 

10.71 

88. 

6.31 

2.37 

2.541 

0 ifi+0.03 
Ui °- 0.03 

19440 

34.41029 

-5.18816 

23.27 

25.17 

10.77 

30. 

4.63 

2.28 

2.291 

n 98+0.02 
0.02 

24828 

34.42530 

-5.16451 

23.30 

24.93 

10.80 

67. 

5.90 

2.27 

2.243 

o 3 1 +0 03 

U-:)1 -0.03 

3956 

34.46006 

-5.25758 

23.50 

26.57 

10.79 

0 

1.46 

2.27 



5326 

34.43100 

-5.25042 

26.08 

27.54 

10.56 

1 

3.19 

2.56 



9277 

34.42780 

-5.23286 

24.41 

26.77 

10.80 

0 

1.61 

2.60 



10771 

34.42257 

-5.22659 

23.55 

25.21 

10.65 

21 

3.73 

2.29 



11700 

34.44719 

-5.22188 

24.57 

26.82 

10.77 

9 

2.71 

2.23 



11909 

34.42215 

-5.22073 

23.88 

27.18 

10.75 

0 

0.63 

2.50 



12447 

34.40580 

-5.21886 

22.79 

25.95 

10.79 

1 

0.86 

2.30 



16478 

34.40424 

-5.20082 

24.33 

27.35 

10.62 

0 

3.97 

2.32 



18367 

34.38511 

-5.19237 

23.60 

27.07 

10.58 

0 

0.76 

2.31 



22984 

34.41577 

-5.17188 

23.81 

26.42 

10.90 

269 

2.44 

2.31 




Notes. 

a 72-band magnitude, defined in Steidel & Hamilton (1992). 
b Sizes calculated in the 77-band as described in van der Wei et al. (2012). 


et al. 2012). The 3D -HST catalogs contain redshifts and stellar 
masses derived from a combination of HST /G141 grism spectra 
and deep photometric data, with wavelength coverage from UV 
to Spitzer / IRAC as described in R. E. Skelton et al. (in prepa- 
ration). Reported stellar masses are the current mass in stars 
and stellar remnants. When possible, the star formation rates 
(SFRs) are based on the UV+1R flux, with the IR determined 
from Spitzer/ MIPS; otherwise, for fainter objects, they come 
from fits of stellar population synthesis models to the stellar 
spectral energy distributions (SEDs). As shown in Wuyts et al. 
(2011) the techniques give consistent answers where they over- 
lap. We note that Ha emission in detected galaxies constitutes 
a median of 5% of the A'-band flux, and thus has a negligi- 
ble effect on broadband photometry. Catalog information and 
emission line properties for the targeted sample are shown in 
Table 1. 

The targets are selected at a fixed cumulative number density 
of 2 x 10 4 Mpc 3 in the stellar mass functions of Marchesini 
et al. (2009), which are ~100% complete in the relevant mass 
and redshift range. Selecting galaxies at a constant number 
density will effectively link galaxies across different redshifts 


if the stellar mass rank order of galaxies is approximately 
conserved with time, or alternatively, if processes that break 
rank order (i.e., merging or scatter in growth rates) have little 
effect on average galaxy properties. Leja et al. (2013) found this 
technique is effective in predicting the median descendant stellar 
mass from z — 3 to z — 0 to ^0.15 dex in a semi-analytical 
model. This selection corresponds to aredshift-dependent stellar 
mass criterion, shown in Figure 1. The stellar mass evolution is 
parameterized as: 


/ M*(z)\ , 

log — LlT =11. 19-0.03z -0.07- 2 . (1) 

\ M q J 


Specifically, we target galaxies that have stellar masses within 
±0.2 dex of this relationship. We further require that Ha and 
[N ii] fall both within the A'-band filter and on the MOSFIRE 
detector, effectively creating a joint constraint on the redshift 
and the (a, S) of each target. These selection criteria result in a 
mask with 16 galaxies, with redshifts 2.1 < z < 2.55. 
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Figure 2. Template SEDs (blue) fit to the broadband photometry (black) are shown for each galaxy with detected emission lines. The units are arbitrarily normalized 
in f\. The center column has //5T/F160W cutouts with the MOSFIRE slit orientation overlaid. The MOSFIRE ^f-band spectra are on the right, with 1 o errors in blue. 
Single Gaussian fits to the Ha and [N n] line profiles are overlaid in red. UDS- 19440 is fit with a double Gaussian to account for the rotational line profile. We also 
include the SED of a typical UV-bright galaxy, constructed with the stellar population synthesis models of Bruzual & Chariot (2003) and typical UV-bright stellar 
population parameters from Papovich et al. (2001). 

(A color version of this figure is available in the online journal.) 


3. OBSERVATIONS AND DATA REDUCTION 

The MOSFIRE A'-band observations were conducted on 2013 
January 20, with ~1"5 seeing. An ABAB dither pattern with 
a 1"5 nod was used. Slit widths were OH. A single mask 
was observed with 16 targets for 85 minutes, with six of the 
targets showing clear line emission. We estimate that sky lines 
obscure only ~3% of the spectral range for emission lines with 
a central per-pixel signal-to noise ratio (S/N) = 20, typical 
for detected emission lines. The non-detections are thus likely 
caused by intrinsically weak or dust-obscured galaxy emission 
lines, rather than overlap of intrinsically bright lines with sky 
emission features. 

The MOSFIRE data reduction pipeline 12 was used to reduce 
the spectroscopic data. The pipeline performs flat fielding, wave- 
length calibration, sky subtraction, and cosmic ray removal be- 
fore producing a final two-dimensional output with an associated 
variance map. One-dimensional spectra were extracted using the 
optimal extraction method of Horne (1986). No flux calibration 
or reddening correction was necessary for this study. 

The A'-band spectra for targets with detected emission lines 
are shown in Figure 2, along with SEDs and F160W direct 
images. Ha and [N n] emission lines are fit with Gaussian 
profiles; the only exception is UDS-19440, which is fit with 
a double Gaussian to properly model the line profile. [N ii] and 
Ha are fit simultaneously, with their line widths and redshifts 
constrained to the same value to improve accuracy when fitting 


12 https://code.google.eom/p/mosfire/ 


the weaker [N ii] line. The adopted fluxes are the areas of the 
Gaussians. Errors in the line profile are determined by perturbing 
each flux value within a Gaussian probability distribution, 
then remeasuring the line profile. The width of the Gaussian 
probability distribution is set to the ler flux error in the pixel. 
The errors on measured parameters are taken as the 68% range 
in derived parameters over 1000 iterations of perturbed spectra. 

4. RESULTS 
4.1. [N n]/Ha Ratios 

We measure a mean [N n]/Ha ratio of 0.27 ± 0.01 (error in 
the mean) ± 0.05 (standard deviation) in our sample. 

We compare the relationship between stellar mass and [N ii]/ 
Ha in our sample and in UV-selected samples (Figure 3). 
By contrasting data rather than derived quantities, we cleanly 
assess potential differences in [N n]/Ha ratio between different 
samples. 

Specifically, we compare our measurements to those of Erb 
et al. (2006a) and Kulas et al. (2013). These studies select 
spectroscopic targets to be bright in the rest-frame UV (1Z < 
25.5) and to fulfill color-color criteria in the rest-frame UV, 
as described in Steidel et al. (2004) and Adelberger et al. 
(2004). The data from Erb et al. (2006a) are stacked spectra 
with ~15 galaxies per point. 

The stellar masses in Erb et al. (2006a) are reported as the 
integral of the SFR. In order to compare with our data, we 
convert the Erb et al. (2006a) stellar masses into the mass in 
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Figure 3. In all panels, blue circles are galaxies with detected emission lines, while black squares are non-detections, (a) Comparing [Nn]/Ha ratios from our 
mass-selected sample (blue) to those from UV-selected samples (orange, red). The mean of our mass-selected sample is marked with a black star. Galaxies with upper 
limits on [N nJ/Hof are marked with downward triangles. The Erb et al. (2006a) data are from composite spectra. The red dashed line is a linear fit to the Erb et al. 
(2006a) data, (b) The UV-bright color-color selection boxes are taken from Steidel et al. (2004). (c) Gray points are from the UV-bright sample of Erb et al. (2006b), 
and the UV-bright magnitude limit is indicated with a dashed line, (d) The size-mass relations are taken from Sersic fits to the circularized H - band light profile at 
Z ~ 2.25 (A. van der Wei, in preparation): median values are shown as solid lines, while the 16th and 84th percentiles are dashed lines, (e) The star forming sequence 
is taken from Whitaker et al. (2012) and marked with a solid red line, with 1 a scatter denoted by dashed lines. 

(A color version of this figure is available in the online journal.) 


stars and stellar remnants using the following formula: 

log ( = 1.06-0.247’ + 0.01T 2 , (2) 

V ^Erb ) 

with T = log (age/yr), and ages taken from Erb et al. (2006a). 
The formula is a fit to the mass loss rates in the Bruzual & 
Chariot (2003) models for a Chabrier IMF, and applies only to 
ages >2 Myr. 


Figure 3 shows that mass-selected galaxies have similar 
[N ii]/Hq: ratios compared to UV-selected galaxies. We quantify 
the significance of this result by simulating our observations 
using population statistics from the UV-selected samples. First, 
we fit a linear relationship to the Erb et al. (2006a) points, 
Ending: 

log ([N ii]/Ha) = -5.36 + 0.44 log (M*/M 0 ). (3) 
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We use this relationship to estimate the UV-selected [Nil]/ 
Ha ratio at the average stellar mass of our sample: log((M}*/ 
Mq) = 10.73. We calculate the biweight scatter (Beers et al. 
1990) about the relationship to be 0.22 dex, using the Kulas 
et al. (2013) galaxies with M* > 10 10 M Q . We then simulate 
our observations by repeatedly sampling six galaxies from a 
Gaussian probability distribution with a mean and scatter fixed 
to the values above. This results in a mean [Nil] /Ha ratio 
greater than our measured ratio only 40% of the time. This 
mass-selected sample thus does not have significantly higher 
[N ii]/Ha ratios than UV-selected samples. 

To further explore whether mass-selected, Ha-detected galax- 
ies are different from UV-selected galaxies of the same mass, 
we include two panels in Figure 3 that explore what fraction of 
the mass-selected sample does not fulfill the selection criteria 
for a UV-bright sample. Our sample is split into galaxies with 
detected emission lines and non-detections. U . Q, and 7 Z mag- 
nitudes are measured directly from the best-fit EAZY template 
(Brammer et al. 2008), then shown relative to the UV-bright 
selection criteria. The results indicate that 50% of the galaxies 
with detected emission lines fit the UV-bright selection crite- 
ria, 1 ' but only 10% of galaxies without detected emission lines 
fit the same criteria. The whole sample is thus primarily UV- 
faint: however, 50% of the galaxies with detected line emission 
are UV-bright. 


4.2. Oxygen Abundances 

We next convert the measured [Nil] /Ha ratios into oxygen 
abundances. To demonstrate the spread in oxygen abundance 
between metallicity calibrations, we calculate the oxygen abun- 
dance in three different calibration systems: Maiolino et al. 
(2008; M08), Denicolo et al. (2002; D02), and Pettini & Pagel 
(2004; PP04). 

The M08 relationship between observed [N ii]/Ha ratio and 
oxygen abundance is: 

log ([N n]/Ha) = co + ci* + C2* 2 + C3* 3 + C4* 4 , (4) 

with co = —0.7732, ci = 1.2357, C2 = —0.2811, C3 = 
-0.7201, c 4 = -0.333, and * = 12 + log(0/H) - 8.69. The 
scatter in this conversion is taken to be 0. 1 dex. 

We find a mean oxygen abundance of 12 + log(0/H) = 
8.87 ± 0.04, 8.70 ± 0. 10 and 8.57 ± 0.08 for the M08, D02, and 
PP04 calibrations respectively. The quoted error is the error in 
the mean, while the standard deviations are 0.09, 0.08, and 0.08, 
again respectively. These mean abundances range from 0.12 dex 
below the solar value of 12 + log(0/H) = 8.69 (Asplund et al. 
2009) for the PP04 calibration to 0. 18 dex above solar for M08. 

4.3. Comparison to Stellar Abundances at z = 0 

We now compare the gas-phase oxygen abundances to the 
stellar oxygen abundances of nearby galaxies. We adopt the 
stellar oxygen abundances measured in Conroy et al. (2013). 
This study analyzes spectra from the inner 0.4-0. 8 effective radii 
of local quiescent galaxies stacked in bins of stellar velocity 
dispersion, and fits a full-spectrum model to them, described 
in Conroy & van Dokkum (2012). The model constrains 
the abundances of individual elements, including oxygen. To 
compare with our results, we derive the average stellar mass 
for these stacks, and interpolate the stellar oxygen abundance 


13 If we instead ask whether these galaxies would be considered UV-bright 
when placed at any redshift, this changes to 67%. 
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Figure 4. Comparing the gas-phase oxygen abundances of galaxies at z ~ 2.3 
to the stellar abundances of local quiescent galaxies at the same number density. 
The mean gas-phase metallicities from multiple metallicity calibrations are 
shown in dark blue. The scatter in these calibrations is ~0. 1-0.2 dex. The solar 
abundance is from Asplund et al. (2009). The gas-phase error bar is the standard 
deviation of the metallicities; it is only shown for the M08 calibration for clarity. 
(A color version of this figure is available in the online journal.) 
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at the expected descendant mass. Specifically, we interpolate 
between the two bins with log(<r /km s^ 1 ) = 2.39 and 2.47, with 
corresponding stellar masses of 10.96 and 1 1.34. This results in 
an oxygen abundance of 12 + log(0/H) = 8.95 ^oq 3 3 . The error 
bars represent the stellar oxygen abundances as inferred at the 
edge of the mass selection box (see Figure 1). 

We can now compare these z = 0 stellar abundances to the 
high-redshift gas-phase metallicities derived in Section 4.2. If 
the high-redshift galaxies are progenitors of the low-redshift 
galaxies, and we are observing the main epoch of star formation 
in the high-redshift galaxies, then the high-redshift gas-phase 
abundance should match the low-redshift stellar metallicity. In- 
terestingly, the gas-phase metallicities are lower than the stellar 
abundances, with the difference depending on the calibration 
method. PP04 produces the largest inconsistency, with the stel- 
lar abundance nearly 0.4 dex higher than the gas metallicity. 
The best match comes from the M08 calibration, which pro- 
duces gas metallicities that are only 0.08 dex lower than the 
stellar abundance. 


5. DISCUSSION 

In this Letter, we have measured [N ii]/Ha ratios from Keck/ 
MOSFIRE K-band spectroscopy of a mass-selected sample at 
high redshift. We demonstrate that these [Nn]/Ha ratios are 
similar to those measured in surveys of UV-bright galaxies. We 
measure an average [Nii]/Ha ratio of 0.27, corresponding to 
an average gas-phase oxygen abundance of 12+log(0/H) = 
8.57-8.87, depending on the metallicity calibration adopted. 
Taking Figure 4 at face value, we would conclude that the M08 
calibration gives the best match to the z = 0 stellar abundances, 
and should therefore be preferred over the other calibrations. 
However, there are many sources of systematic uncertainty in 
this study. 

We first consider the sources of the differences in the three 
calibrations. The M08 oxygen abundance is calculated from a 
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combination of the direct T e method for galaxies with log(0/H) 
+ 12 < 8.35, and with strong line ratios and the photoionization 
models of Kewley & Dopita (2002) for log(0/H) + 12 > 8.35. 
D02 derives oxygen abundances primarily from the direct T e 
method, with roughly one third of their abundances from oxygen 
or sulfur strong line ratios. PP04 uses 97% direct T e abundances. 
All three studies then derive a relationship between [Nil] /Ha 
ratio and oxygen abundance. 

The spread in metallicity calibrations may be related to issues 
with the photoionization models, or temperature gradients and 
fluctuations that cause metallicities determined by the direct 
T e method to be underestimated (Peimbert 1967; Kewley & 
Ellison 2008). Observational studies of optical recombination 
lines show a systematic differences of 0.26 ± 0.09 dex with the 
T e method (Esteban et al. 2009), which may explain some of the 
variation. 

There is also emerging evidence that, at high redshift, metal- 
licities based on [N n]/H« are offset from those based on 
[Oii]A3727, [Oiii]A4959, 5007 and H/l, even within the same 
system of metallicity calibrations (Newman et al. 2013; Cullen 
et al. 2013). Specifically, metallicities based on oxygen lines are 
~0.3 dex lower than those based on nitrogen lines. This differ- 
ence is attributed to different physical conditions in star forming 
regions at high redshift, and it is unclear which calibrations, if 
either, represents the “true” gas-phase metallicity. 

Another uncertainty in our study stems from the fact that, in 
addition to metallicity, [N ii] flux is sensitive to the presence of 
both an active galactic nucleus (AGN) and shock excitation. 
It is unfortunately not possible to separate out the AGN 
contribution to the [N ii] flux without high-resolution integral 
field unit (IFU) data (Forster Schreiber et al. 2011). However, all 
detected galaxies in our sample have log([N n]/Ha) < —0.5; the 
theoretical high-redshift Baldwin, Phillips, & Terlevich diagram 
indicates that significant AGN contribution in this regime is 
unlikely (Kewley et al. 2013), though it may exist at a low level. 

We note that the presence of radiative shocks, AGN, and/or 
different physical conditions at high redshift would imply lower 
oxygen abundances than indicated by the measured [Nil] /Ha 
ratios. 

It may also be that the stellar oxygen abundances are over- 
estimated. The measurement of oxygen abundance in old un- 
resolved stellar populations is notoriously difficult (see Conroy 
2013 and references therein). In the Conroy & van Dokkum 
(2012) spectral model, the oxygen abundance is primarily de- 
rived from a combination of TiO lines and molecular equilib- 
rium involving CNO. The trend [0/Fe]~[Mg/Fe] matches that 
as measured in the Milky Way (Edvardsson et al. 1993, though 
see also Bensby et al. 2013), suggesting that the stellar oxygen 
abundance measurements are robust. 

We also consider the possibility of misidentification of the 
descendant galaxies. The predicted stellar mass evolution from 
z ~ 2.3 to z = 0 for these galaxies is +0.4 dex. If instead 
the stellar mass evolution were +0.1/+0.7 dex, this results in 
negligible change of — 0.05/+0.05 dex in the oxygen abundance 
of their descendants. This is due to the flat relationship between 
the stellar mass and stellar metallicity (Panter et al. 2008). This 
comparison is thus robust against errors in descendant matching. 
Additionally, the stellar mass growth at this number density 
is primarily not through additional star formation, but rather 
accretion (van Dokkum et al. 2010; Patel et al. 2013). These 
accreted galaxies presumably had lower metallicities than the 
main progenitor galaxy, which means that they would lower 
the total luminosity-weighted stellar metallicity (Greene et al. 
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2013). Taking this process into account would therefore, if 
anything, increase the discrepancy in Figure 4. 

Perhaps the most likely possibility is that selection by stellar 
mass may still be biased, in the sense that the most metal-rich 
galaxies do not have detectable emission lines. Only 37.5% 
of our sample has detected emission lines; furthermore, the 
galaxies without detected emission lines are not a random 
subset of the sample. The last two panels in Figure 3 examine 
the different properties of detections and non-detections in our 
sample. Galaxies with detected emission lines consistently have 
larger effective radii and higher SFRs than galaxies without 
detected emission lines. 

These properties are known to correlate with the gas-phase 
metallicity of the galaxies. For example, observations suggest 
that gas-phase metallicity and SFR are anticorrelated (Mannucci 
et al. 2010). Furthermore, gas-phase metallicity varies by up to 
0.2 dex at fixed stellar mass as a function of half-light radius, 
with larger galaxies having lower metallicities (Tremonti et al. 
2004; Ellison et al. 2008). 

Even in galaxies with detected emission lines, some unknown 
fraction of the star formation may be obscured by dust. Perhaps 
the detected line flux originates from “shells” or “rings” of 
Ha, as seen in IFU studies (Forster Schreiber et al. 2011), 
while the star forming core remains heavily obscured. UDS- 
22984 may host such obscured star formation: it has a 24 /im 
flux indicating a SFR of ~270 M 0 yr” 1 and an SED-estimated 
A v = 2.6, yet no detected line emission. Since dust correlates 
with metallicity, this galaxy will likely be more oxygen-rich 
than galaxies with unobscured star formation. If such heavily 
obscured star formation is common at high redshift, it remains 
a possibility that the star forming progenitors of local ellipticals 
have yet to be detected. 
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